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Abstract 

The NA60 experiment at the CERN SPS has studied dimuon production in 158 AGeV In-In 
collisions. The strong pair excess above the known sources found in the mass region 0.2 < 
M < 2.5 GeV has been previously interpreted as thermal radiation. In this paper results on 
the associated angular distributions for M < I GeV, as measured in the Collins-Soper reference 
frame, are presented. The structure function parameters A, fi, v are consistent with zero and the 
projected polar and azimuth angle distributions are uniform. The absence of any polarization is 
consistent with the interpretation of the excess dimuons as thermal radiation from a randomized 
system. 

Lepton pairs are a particularly attractive observable in high energy nuclear collisions, because 
their continuous emission probes the entire space-time evolution of the produced fireball. To the 
extent that the bulk consituents of the expanding matter (hadrons and partons) equilibrate, the 
directly generated lepton pairs should appear as a "thermal radiation". Such a thermal radiation 
should exhibit a number of features: (i) a Planck-like exponential shape of mass spectra (strictly 
correct in case of a flat spectral function), (ii) exponential mr spectra, (iii) the absence of any 
polarization in the angular distributions. 

Previous work on the mass spectra showed that dilepton production for M < 1 GeV is largely 
mediated by the process kk p —> /u/u with a strongly broadened p (7]. The analysis of 
the pt spectra shows an increase of T^ff vs M up to ~ 1 GeV, while the drop at ~ 1 GeV 
and above might signal the transition to the partonic thermal process qq — > y //// |0 IH. 
The acceptance-corrected mass spectrum falls exponentially and it is described quantitatively by 
theoretical models which are explicitely based on medium thermalization up to 2.5 GeV [3]. 

This paper concentrates on the first measurement of dilepton angular distributions in high 
energy nuclear collisions. Results were already recently published in [4]. Here, further details 
and aspects on the results will be provided. In general, the diff'erential decay angular distribution 
in the rest frame of the virtual photon with respect to a certain set of axes can be written as Isj] 



1 d^o- 



cr dcos 



-^l + Acos^ 6 + jj sin 20 cos + ^ sin^ 0cos20j . (1) 



Here A, /i and v are structure functions related to helicity structure functions and the spin 
density matrix elements of the virtual photon. The nomenclature follows 0]. Because of 
lack of sufiicient statistics for M > I GeV, the present analysis is restricted to M < 1 GeV, 
which is dominated by pion annihilation. Even in that case, where the annihilating particles 
are spinless, the structure functions can be different from zero because of orbital momentum. 
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However, a completely random orientation of annihilating pions in 3 dimensions, as expected in 
a thermalized medium, should lead to A, ij,v - 0. 

The decay angular distributions were studied in the ColUns-Soper (CS) reference frame J^. 
The z axis is the bisector between the beam and the negative target momenta (defining the reac- 
tion plane). The polar angle 6 is the angle between the positive muon and the z-axis (defining 
the decay plane), while (p is the angle between the reaction and decay planes. The choice of the 
frame is not relevant: once all measured. A, n and v can be re-computed in any other frame with 
a simple transformation fi"]. This point will be illustrated in the results. 

Details of the NA60 apparatus are contained in fllst]. The data sample for 158 AGeV col- 
lisions is the same as in |1, 2]. The analysis was done in the 2-dimensional cos0 - space 
with different binnings in (d^N/d cos 0d<p)ij in order to assess the systematics. For each i,j bin, 
the analysis proceeded with the following steps: (i) subtraction of the combinatorial background 
from n and K decays by the event-mixing technique, (ii) subtraction of fake matches (incorrect 
associations of muon tracks to the tracks reconstructed in the silicon pixel vertex telescope), (iii) 
isolation of the excess by subtraction of the known sources, (iv) acceptance correction in the 
2-dimensional bin. Further details for the analysis steps can be found in IH |3, S 0] • In order 
to increase the statistical significance, no centraUty selection was required. In order to exclude 
the region of the low mj rise seen at all masses |12|], a cut pr > 0.6 GeV was applied. That has 
also the effect of improving the background/signal ratio, which is ~2-3. The excess dileptons 
were studied in two mass windows: 0.4 < M < 0.6 GeV (~ 17600 /i/i pairs) and 0.6 < M < 0.9 
GeV (~ 36000 yUyU pairs). A fully differential acceptance correction should be performed in the 
5-dimensional space M - pj -y - cos 0-(p. In this analysis the Monte Carlo was first tuned with 
an iterative procedure to the measured M, y and pj spectra and then the acceptance matrix in the 
cos 9 - (p space was determined. The Monte Carlo muons were overlayed on real data to include 
the effects of pair reconstruction efficiencies. 




Three different methods were applied to extract the structure functions. The 2-dimensional 
distribution d^N/d cos 9d(p was directly fitted with Eq. (1). In order to avoid very low acceptance 
bins, the 2-dimensional distribution was restricted to a 6x6 matrix with -0.6 < cos < 0.6 and 
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-0.75 < cos4> < 0.75. Alternatively, fixing yu to the value found in the 2-dimensional fit, the 
2-dimensional distribution can be projected in cos 6 and (p and these 1 -dimensional projections 
can be fitted with dN/d\ cos9\ oc (l + Acos^ 6^ and dN/d\(p\ oc (1 + A/3+v/3 cos 20), respectively 
(method 2). Finally, an analysis of the inclusive distributions in | cos 6\ and \<p\ can be also per- 
formed (method 3). In the analysis of 1 -dimensional projections (methods 2 and 3), the range 
in and </> can be enlarged to | cos 6\ < 0.8, < | cos <p\ < 0.75 for method 2 and | cos 6*1 < 0.8, 
< 101 < TT for method 3. 
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Figure 2: Polar angle distributions for excess dileptons. 
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Figure 3: Azimuth angle distributions for excess dileptons. 



The left panel of Fig. [T| shows the 2-dimensional contour plot from the fit for the excess 
dileptons in the two mass intervals 0.6 < M < 0.9 GeV, while the right panel shows the cor- 
responding projections in cos 6 for the different bins in cp. Similar results are obtained for the 
mass interval 0.4 < M < 0.6 GeV. The parameter fi, which can be determined only from these 
bidimensional fits, is consistent with zero in both cases: yu = 0.05 + 0.03 for 0.6 < M < 0.9 GeV 
and fi - -0.04 + 0.10 for 0.4 < M < 0.6 GeV. The parameters A and v are determined with a 
better statistical accuracy from the analysis of the 1 -dimensional projections where the low pop- 
ulated bins in the acceptance borders have a smaller effect. The polar and azimuth distributions 
obtained from methods 2 and 3 are shown in Fig. |2] and Fig. |3] respectively. Also A and v are 
consistent with zero for both mass windows. 

As already mentioned, once all measured. A, jj. and v can be re-computed in any other frame. 
Fig.|4]shows the values of A and v determined from the fits in the CS frame and in the Gottfried- 
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Figure 4: Structure functions A and v in tlie Collins-Soper and Gottfried- Jackson frames. 

Jackson frame as obtained by applying the transformation of [6]. The structure functions remain 
consistent with zero independent of the frame. 

The systematic error due to the combinatorial background subtraction is ~ 2 - 3%. The 
systematic error due the fake matches is less than 1%. Variations in M, y and pr used as input 
for the Monte Carlo lead to uncertainties much smaller than the statistical errors. Systematic 
uncertainties in the isolation of the excess range from 4 - 6% up to 10 - 15% in some low- 
populated cos 6 - (f) bins and they represent the main source of systematic errors. However, 
the measurement still remains dominated by the statistical errors. If the systematic errors are 
conservatively assumed uncorrected from point to point and added in quadrature to the statistical 
fit errors quoted in the numbers above, this would increase by 15-20%. 

Summarising, the absence of any polarization is fully consistent with the interpretation of 
the observed excess as thermal radiation. It is important to remember that this is a necessary 
but not sufficient condition. Put together with the other features - Planck-like shape of mass 
spectra, temperature systematics, agreement of data with thermal models - this makes the thermal 
interpretation more plausible than ever before. 
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